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Abstract 

We extend the quark Born diagram formalism for hadron-hadron scattering to pro- 
cesses with valence qq annihilation, specifically 1 = and 1 = 1 vrvr and / = 1/2 
Ktt elastic scattering. This involves the s-channel hybrid annihilation process q^q^ — > 
qqg -^ q^(f and conventional s-channel qq resonances, in addition to the t-channel 
gluon exchange treated previously using this formalism. The strength of the t-channel 
gluon amplitude is fixed by previous studies of / = 2 vrvr and / = 3/2 Ktt scattering. 
The s-channel resonances /9(770), K*{S92), /o(1400) and Ko(1430) are incorporated as 
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relativized Breit-Wigner amplitudes, with masses and energy-dependent widths fitted 
to experimental phase shifts. The strength of the s-channel gluon "hybrid" annihilation 
diagrams is problematical since the perturbative massless-gluon energy denominator 
must be modified to account for the effect of confinement on the energy of the vir- 
tual hybrid state. Our naive expectation is that near threshold the hybrid diagrams 
are comparable in magnitude but opposite in sign to the contribution predicted using 
massless perturbative gluons. Fitting the strength of this amplitude to / = vrvr and 
/ = 1/2 Ktt S-wave data gives a result consistent with this expectation. We find good 
agreement with experimental phase shifts from threshold to 0.9 GeV in tttt and 1.6 
GeV in Ktt using this approach. We conclude that the most important contribution to 
low energy S-wave scattering in these channels arises from the nonresonant quark Born 
diagrams, but that the low-energy wings of the broad s-channel resonances /o(1400) 
and i^Q(1430) also give important contributions near threshold. The nonresonant con- 
tributions are found to be much smaller in L > partial waves, but may be observable 
in / = 1 TTTT and / = 1/2 Kn P-wave phase shifts. 



1. Introduction 

The determination of hadron-hadron interactions in terms of quark and gluon degrees of 
freedom is an important goal of the nonrelativistic quark potential model. Historically these 
studies first concentrated on the nucleon-nucleon interaction, which is well determined ex- 
perimentally and is fundamental to nuclear physics. Nonperturbative techniques such as the 
resonating group method [jl| and variational approaches lead to reasonable descriptions 
of the short- and intermediate-range nucleon-nucleon interaction in terms of quark forces. 
(For a review of this work on NN to 1989 see Shimizu 0.) Studies of other hadron-hadron 
channels in terms of quark and gluon interactions have also been successful. Weinstein 
and Isgur 0] carried out variational studies of pseudoscalar meson-meson interactions us- 
ing the nonrelativistic quark potential model, and found weakly-bound deuteron-like ^^KK 
molecules", which have been identified with the /o(975) and ao(980). They also used their 
variational techniques to extract equivalent low-energy meson-meson potentials, which with 
some qualifications give reasonable results for tttt and Kn elastic phase shifts [^, |^. The 
resonating group techniques have also been applied to the kaon-nucleon system and give 
results similar to the observed 1 = and 1 = 1 S-wave phase shifts. There may be some dis- 
crepancies, however, and the 1 = KN phase shift in particular is not yet well determined 
experimentally and merits further study, which may be possible at DAPHNE 0. 

Recently Barnes and Swanson 0, |l^ introduced a perturbative Born-order formalism 
for hadron-hadron scattering in terms of quark and gluon degrees of freedom, based on a 
single interaction, usually one-gluon-exchange (OGE) followed by quark line interchange. 
This analytical technique reproduces many of the successes of the earlier nonperturbative 
calculations when applied to scattering processes which are free of valence quark annihilation. 
In some cases such as the / = 2 tttt and / = 3/2 Kn S-waves the Born diagrams are in 
remarkably good agreement with experiment from threshold to the maximum experimental 
invariant mass of about 1.5 GeV. The reactions studied to date using this method are 1 = 2 
TTTT [|, |Tg], / = 3/2 Ktt [|n[, certain vector- vector channels [|TU|, 0, / = 0,1 KN 



and NN, NA and A A |T^, and several other related channels. One conclusion of these 
studies was that the powerful but complicated nonperturbative techniques were unnecessary 
in some channels, notably vrvr and Kn, because the perturbative amplitude alone gives a 
good description of the data. 

Motivated by the successes of this simple perturbative approach, in this paper we gen- 
eralize this technique to scattering processes with valence annihilation. General hadron- 
hadron scattering amplitudes include contributions from resonance production, resonance 
exchange and qq annihilation (s-channel gluon exchange), in addition to the t-channel gluon 
exchange considered previously in this formalism. In principle each of these contributions 
may be important, so reactions in annihilation channels can be quite complicated. Here we 
consider pseudoscalar-pseudoscalar (Ps-Ps) elastic scattering with valence annihilation, and 
assume that the important amplitudes arise from t-channel gluon exchange, s-channel gluon 
exchange and s-channel resonance production. The first contribution is treated using the 
quark Born diagram techniques developed previously, and we find that the t-channel gluons 
which dominated J = 2 tttt and / = 3/2 Kit are numerically rather weak here. The sec- 
ond scattering mechanism, s-channel gluon exchange, is treated using the same quark Born 
techniques. Special care must be taken with this process, however, to incorporate the mass 



of the intermediate hybrid state. This changes s-channel gluon exchange from a repulsive 
to an attractive interaction in the cases studied here, and makes the overall magnitude of 
this effect rather uncertain. Here we fit it to experimental S-wave phase shifts, and compare 
the fitted strength to our naive estimate based on an effective gluon mass of ~ 1 GeV. The 
third contribution, from s-channel resonances, is incorporated phenomenologically by treat- 
ing these as relativized Breit-Wigner phase shifts with free masses and widths, which we fit 
to the data. 

In section 2 we describe our techniques in detail for / = and / = 1 tttt scattering, and 
give numerical estimates of the contribution of each effect to the scattering length. In section 
3 we give the corresponding results for J = 1/2 Kti scattering, which is formally similar to 
TTTT but is somewhat more complicated due to the strange quark mass. In section 4 we carry 
out detailed fits to experimental data sets for J = and / = 1 tttt and 7 = 1/2 Kti phase 
shifts and discuss the relative importance of resonant and nonresonant contributions. We 
conclude with a brief summary of our results and suggestions for future work. 

2. Detailed Results for 7 = and / = 1 tttt Scattering Amplitudes 

We consider three contributions to hadron-hadron scattering in annihilation channels, 

1) t- channel gluon exchange, 

2) s-channel gluon exchange (valence qq annihilation), and 

3) s-channel resonances. 

A fourth contribution which is often discussed is t-channel meson exchange. We exclude this 
scattering mechanism in tttt and Ktc scattering because the most important such contribu- 
tion, one pion exchange, is not present in Ps-Ps scattering, and in any case the contribution 
of t-channel meson exchange to low-energy hadron-hadron scattering has probably been 



overestimated. For a discussion of this issue see |]T5 . 

Taking the three mechanisms in order, their contributions to the scattering amplitude 
are as follows: 

1) t-channel gluon exchange 

In earlier studies of Ps-Ps elastic scattering [4-6] it was found that the dominant scattering 
mechanism involves the OGE spin-spin hyperfine interaction, which between quarks i and j 
is 

Sniimj 

A similar conclusion was reached earlier for NN scattering [1-3]. An explanation of hyperfine 
dominance in Ps-Ps scattering was given by Barnes and Swanson |P, |lOl, who found that 
the matrix elements of A" ■ X'j Si ■ Sj in the four Born-order quark scattering diagrams for 
{qq){qq) — > {qq){qq) all have the same sign in this channel. In contrast, the spin-independent 
matrix elements of A" ■ A" (which multiply the spin-independent color Coulomb and linear 
confining terms) interfere destructively between the diagrams. In the Ps-Ps annihilation 



H^, = -Tj—^ (A- /2) ■ (A,72) {S, ■ S,) 5{n,) . (1) 



channels we consider here the same conclusions apply, so we again assume dominance of the 
spin-spin color hyperfine term in t-channel gluon exchange. 

The t-channel gluon exchange contribution to J = vrvr scattering can be determined 
from the previous quark Born study of / = 2 tttt scattering ||^. The only difference is the 
flavor factor, which gives a Hamiltonian matrix element of 

h^fT^{t—ch.. gluon) = — -/ij^^(t— ch. gluon) = 



Alias 1 
9m2 (27r)3 
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The J = 1 vrvr matrix element from t-channel gluon exchange is zero; this result and the 
/ = 0// = 2 ratio both follow immediately from the observation that vr+vr^ -U vr+vr^ through 
this mechanism. Here /z = cos(^c.m.), where 6c,m. is the center of the mass scattering angle, 
and k is the magnitude of the asymptotic three- momentum of each meson in the cm. frame. 
This is the matrix element of the spin-spin OGE term calculated between single- Gaussian 
qq wavef unctions, summed over all four permutations of gluon exchanges and quark line 
interchanges. The wavefunction parameter jSt^ is usually taken to be about 0.3 GeV in the 
nonrelativistic quark model. In a fit to the / = 2 vrvr S-wave phase shift gave (3.,^ = 0.337 
GeV, as = 0.6 and niq = 0.33 GeV, which we also use here. (Note that the scattering 
amplitude (2) from the spin-spin interaction only involves the combination as/Tn"^.) 

The Hamiltonian matrix element hfi is proportional to the Born-order T-matrix element 
and can be used to calculate Born-order phase shifts through the relation 



^bIvu = ^^ J_^ hfi{n)Pe{fi)dfi , (3) 

where Ei and E2 are the two hadron energies in the cm. frame (set equal for vrvr), E is the 
total cm. energy Ei + E2, and iS is a statistical factor which is 2 for tttt and 1 for Kit. The 
resulting / = vrvr phase shifts for even i are 



4orn """^"(i-ch. gluon) = +^kE^ 



V27 



(4) 



where x = /c^/4/5^. The S-wave phase shift, using io{x) =sinh(x)/x, leads to a scattering 
length of 



ar (t-ch. gluon) = hmf^yk = +9 (,1 + yf^j ^ • (5) 

This phase shift and scattering length are positive, corresponding to an attractive inter- 
action, but are numerically rather small; with our parameters this scattering length is 
aQ^° = -1-0.043 fm, an order of magnitude smaller than the experimental value of 

ni=0(^^r.i ^ - / +0-32(13) fm (production expts), , . 

«o ^^^P^-^ - \ +0.37(7) fm (ire4 decay). ^^' 

These numbers are taken from a recent review by Ochs [0, and incorporate constraints 
from dispersion relations. 



This conclusion regarding the small contribution of non-annihilation processes in / = 



TT7T is more general than this model, since a^ 



2^0 



follows from tt+tt -/^ tt~^7; in any 



single-pair-interchange model, and one can use the experimental / 

of Cg^^ ~ —0.08 fm [Il7| to normalize the no-annihilation amplitude. 



2 TTTT scattering length 
Evidently S-wave 



J = vrvr scattering is dominated by annihilation processes, and the interesting question is 
whether these are due primarily to a broad /o qq state (as is often assumed) or to a broad 
scalar glueball or another intermediate state such as qqg. 

2) s-channel gluon exchange 



The s-channel gluon exchange scattering mechanism is in many ways the most interesting. 
We will first calculate this contribution using perturbative gluons. If gluons behaved like 
photons we could model the low-energy effects of s-channel gluon exchange by the familiar 
positronium annihilation interaction W^ augmented by a color factor. 



Hf''\s-di. gluon) = 4 



27ra, 



m„ 



(A?/2) ■ (A^/2) (5. ■ S, + -) 5{n 



(7) 



For s-channel transverse gluon exchange this is the leading term in an expansion in v"^ /(?. 
There is also a color Coulomb interaction, but we expect this to be small because the 
annihilation color charge density (0|/9"(x)|gg) transforms as L = 1. 

A typical annihilation diagram involving this Hamiltonian, for K^-k^ -^ K^tt' , is shown 
in Fig.l. The contribution of this diagram to the meson-meson scattering amplitude can be 
derived easily using the diagrammatic techniques discussed in , see especially Appendix C 
of that reference. In Ps-Ps scattering this Hamiltonian has a color matrix element of +4/9 
and a spin matrix element of +3/4. We may separate the isospin amplitudes by considering 



the special cases tc^tc -^ tt+tt and tt+tt" 



n TX ' 



(^/( 



TT- 



— ^/(Tr )); this leads to 



/ijj °(s— ch. gluon) = 4 



27ras 1 
m^ (27r)' 



and 



/.-^(^-ch. gluon) = +^ ^ 



g-(i-M)a; _|_ ^-{l+^l)x 



^-{l-^l)x _ ^-{l+^l)x 



(9) 



again using x = k'^/A/S^. The 1 = 2 annihilation amplitude is of course zero. Using (3), these 
amplitudes lead to / = and / = 1 tttt phase shifts of 



<^Born('S-ch. gluou) = C/,£ -^ kE^ 6 ^ if,{x 



mt 



(10) 



where C/^£ is —1/2 for / = and i =even, —1/3 for /= 1 and 
The S-wave phase shift gives an / = scattering length of 

7=0/ , 1 X IttsJ^TT 

% (s-ch. gluon) = — . 

2 mi 



=odd, and zero otherwise. 
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This negative phase shift and scattering length correspond to a repulsive interaction, opposite 
to the observed / = vrvr interaction. With our standard quark model parameters this s- 
channel massless-gluon contribution is Oq^^ = —0.076 fm, somewhat larger than t-channel 
gluon exchange and opposite in sign. 

More careful consideration suggests that the assumption of massless perturbative gluons 



is especially unphysical in this case and should be modified [|l9l . One effect of confinement 
is to raise the mass of the lowest intermediate qqg hybrid basis state from the perturbative 
value of 2mg to the physical hybrid mass M^, which is well above the vrvr invariant masses we 
consider. This causes the energy denominator 1/ {Egg — Eg) in the qq -^ g ^ qq subprocess to 
change sign. Since the sign of the amplitude for this second-order process is also determined 
by the energy denominator, confinement changes this sign as well, which leads to an attractive 
force in / = TTTT and / = 1/2 Kit. There are many ways one might incorporate this effect 
of confinement on virtual gluons; one simple approach is to modify the denominator of the 
gluon propagator by including an effective gluon mass, 

s-'^{s-fiiy\ (12) 

where the mass scale fig is set by the hybrid mass gap, which is presumably about 1 GeV. 
This modification of the gluon propagator changes the effective low-energy qq interaction 
in (7) by replacing ttt,"^ by [m'^ — /iV4)^^. This multiplies the naive Hamiltonian (7) by a 
negative constant, 

i/f ""'"*■ (s-ch. gluon) = -I - — /2^ \2 _ X f ' Hr'^^'is-ch. gluon) . (13) 

If we estimate constituent masses of fig = 1. GeV and nig = 0.33 GeV, we expect the s- 
channel effective interaction to be similar in magnitude to the original massless gluon form 
but opposite in sign. The exact numerical strength of this diagram is clearly problemati- 
cal due to uncertainties in the effect of confinement, so in practice we simply introduce a 
multiplicative factor / in the s-channel annihilation Hamiltonian and use the form 

Hr^\s-ch. gluon) = +/ ■ ^ (A?/2) ■ (A^/2) (5, • 5, + h 6(f,,) . (14) 

nig 4 

The phase shifts (10) and scattering length (11) from s-channel gluon exchange are thus 
multiplied by a phenomenological parameter /. We will fit / to the experimental J = tttt 
and 1 = 1/2 Ktt S-wave phase shifts, with the caveat that we expect it to be negative and 
not far from unity. 

3) s-channel resonances 

We treat the phase shift due to conventional s-channel qq resonances phenomenologi- 
cally, since it is not clear how to determine three-meson couplings directly from quark-gluon 
interactions. We use a relativized Breit-Wigner form suggested by the Particle Data Group 
20| to model the s-channel resonances, which has an elastic phase shift of 



^— if^}. <-) 



where s = E"^. Since we are considering broad resonances it is important to include the energy 
dependence of the width T{E). Of course this function is somewhat model dependent. In 
general we expect it to consist of a centrifical factor for decay to a two-body final state with 
orbital angular momentum L, times a phase space factor of k/M and a vertex form factor 
D{k) derived from wavefunction overlaps of the initial and final mesons. Here we use the 
form suggested by the LASS collaboration [^], 

In this formula kfi is the momentum of one final-state hadron at the s-channel resonance 



mass, and for vrvr, k^ = J Mjj^/A — m'^. (Here and in the remainder of the paper we assume 
single decay channels for each resonance, tttt for /q and p and Kit for K^ and K*. The 
broad scalar resonances couple dominantly to these channels so this is a reasonable first 
approximation, especially for low-energy effects.) 

The vertex form factor D{k) is wavefunction- dependent; here we use 

D{k)=cexp\-kyQf3l\ , (17) 



which follows from the ^Pq decay model with SHO wavefunctions |2^. We expect the 
corresponding form factor for radial excitations to suppress the low-energy contributions 
of radially-excited qq states so they are not important numerically. 

The low-energy limit of this resonance phase shift gives the contribution of a scalar 
resonance to the scattering length, which in the vrvr (equal mass) case is 

ao(res. 



(1 - 4m2 /M2)3/2 Dikn) Mr Mr' 

For 7 = TTTT scattering the only well established broad resonance is the /o(1400). The 
mass and width of this broad state are both rather problematical; the PDG mass and width 
estimates are Mr ^ 1400 MeV and T{Mr) ^ 150 — 400 MeV. For this wide range of values 
the scattering lengths with the ^Pq form factor are 

ao(res.) = +(0.06 -0.17) fm . (19) 

Taking the recent Crystal Ball central values for the mass and width from 77 -^ 7r°7r°, M = 
1250 MeV (assumed) and r_R= 268(70) MeV [^, we find a contribution to the scattering 



length of 

ao(res.) = +0.11(3) fm . (20) 



The Crystal Barrel collaboration |£4[ recently reported a similar mass and width for this 
broad /o, M ^ 1335 MeV and F = 255(40) MeV, corresponding to a very similar scattering 
length of 0.11(2) fm. The Crystal Barrel and Obelix [^ collaborations report a similar, 
somewhat broader /o in An final states, with M = 1374(38) MeV and F = 375(61) MeV 
and M = 1345(12) MeV and F = 398(26) MeV respectively. These imply a resonance 
contribution to the scattering length near the upper limit of (19). 



Recall for comparison that the nonresonant gluon exchanges give contributions to the 
scattering length of 

ao(gluon ex.) = { + 0.043 - 0.076/ } fm , (21) 

where / is expected to be negative and not large relative to unity. Since the experimental 
scattering length is about +(0.3 — 0.4) fm (6), these results suggest that low energy S-wave 
/ = TTTT scattering receives important contributions both from the broad /o(1400) and from 
nonresonant scattering, and that / ^ —2 to —3. Our detailed fits to phase shifts will support 
similar values for /. 

The resonance contributions remain large at low energies because of the 1/ E factor in 
the energy-dependent width T{E) (16). Although this appears well motivated as the 1/M 
in the decay rate of an initial state of mass M [^, ^], this factor and the (5(r, E) relation 
(15) are so important to the threshold behavior that they merit more careful study. It would 
also be interesting to explore the sensitivity of our conclusions to the form chosen for the k- 
dependent three-meson vertex; using instead a point like-meson form factor, we would predict 
a scattering length from the broad /o resonance about half as large. This is because the form 
factor suppresses the coupling of the /q to pions at higher momenta, so if we use the width 
at resonance as a fixed input, the strength of the coupling at threshold is increased by the 
form factor. (Note the D{0)/D{kji) dependence of the scattering length in (18).) Another 
concern is that scattering amplitudes from two different time-orderings, vrvr — > /o -^ vrvr and 
vrvr — i> /ovrTTTrvr -^ tttt, are added with equal strength to give the covariant form (15). These 
are actually modified by form factors, and the second "Z-graph" process may be strongly 
suppressed p8| . We will test the importance of some of these complications in our detailed 



study of experimental S-wave phase shifts. 

3. / = 1/2 Kit Scattering Amplitudes 

Application of these techniques to / = 1/2 Kit scattering is straightforward. First, t-channel 
gluon exchange for I = 1/2 is simply related to the / = 3/2 amplitude, 

h'='^\t-ch. gluon) = -i h'jf/'it-ch. gluon) , (22) 

which follows from K^tt^ -/-^ K^n^ through this mechanism. Taking the / = 3/2 result 
from [|TI|, we have 



h'='^\t-ch. gluon) = - (^^ (Ti +T, + C, + C) , (23) 

where Tj and C, represent the contributions of the "transfer" and "capture" diagrams 0. 
Specializing to the case of identical pion and kaon spatial wavefunctions {(3.„ = Pk), these 
are explicitly 

T.^exp(-i^±^^±pil^4, (24) 



2-2C + C' + 2(l-C)/i , .^,. 

T2 = p exp <^ X >, (25) 



Ci = p (4/3)"'/ exp| x\ (26) 

and 

C. = (4/3Pexp{-i±i±^i!^(5i±^.}. (27) 

where x = k'^/4j3'^, as in vrvr scattering. These somewhat comphcated results allow for 
different strange and nonstrange quark masses through the parameter 

p = nig/ms , (28) 

and ( is the combination (1 — p)/(l + p). In our previous study of / = 3/2 Ktt scattering 
pj| we considered the more general problem of different pion and kaon length scales, but 
found that the I = 3/2 scattering amplitudes were quite insensitive to the ratio Pn/Px, and 
setting it equal to unity allowed a very good description of the data. For equal length scales 
we found an optimum value of p = 0.677, which we will also assume here. 

The s-channel gluon exchange diagram with /5^ = Pk gives a Hamiltonian matrix element 
ftr'^(-ch.gluon) = +/.^J:^exp(-(i±ii±i!lii^4. (29) 



''fi V' ^^^- &iuuii; — -rj 2i^ 

' m^ [ZiT 

which by inspection is proportional to the contribution of the t-channel gluon diagram 
"transferi" (Ti above), and so can be incorporated in (23) by the substitution Ti -^ (1 — 
9//2)Ti. The general-^ phase shift and the corresponding scattering length from combined 
s- and t-channel gluon exchange are 



^1=1/21/ 1 X o^s kE^EK 

' ' (gluon ex. 



9ml E 



9^N ,^{i+(+he)x : /2 + 2C + C' 



_;)e-(i+c+tc-)^,^^ - - - X 



■i(4+c+2c'')^ ,• r^^jLL 



+p {e-(2-C+K^)-+ (4/3)^^/2 g-i(4-c+C^). I ,^(^a;) + (4/3)3/2 g-i(4+c+2C^)- ^,(^^r^ x 
and 



7=1/2/ , N tts m^rriK 

Oq (gluon ex.) — 



(l + (4/3)^/^)(l + p)-^/ 



(30) 



(31) 



9ml m,T^ + wlk 

The scalar iirQ(1430) is the single important s-channel resonance in low-energy S-wave 
1 = 1/2 Ktt scattering. We incorporate this state using a relativized Breit-Wigner of the 
form (15-17), as in vrvr scattering. It leads to a scattering length of 

«o ires.j ^2_(^^ + ^^)2 k^ ^(^^) Mn Mn ^ ^ 

The PDG mass and width for this state (taken from LASS results) are M = 1429±4±5 MeV 
and r = 287 ± 10 ±21 MeV [^], but subsequent reanalysis by this group has led to preferred 
values of M = 1412 MeV and F = 294 MeV [^]. These correspond to a contribution to the 
scattering length of ao(-ft'o(1430)) = 0.145 fm given our form for the relativized Breit Wigner 
(presumably with a statistical error from F of about 10%). In summary, the total scattering 
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length we expect from gluon exchange and the s-channel resonance i^Q(1430) is numerically 
(excluding errors) 

aJT^^^iKTT thy.) = j 0.056 - 0.059/ + 0.145| fm . (33) 

which we may compare with the experimental scattering length found by Estabrooks et al. 



3g 



ai'^^^Kn expt.) = 0.472(8) fm . (34) 

Evidently the scattering lengths (33) and (34) suggest a value of / ~ —4 for the s-channel 
gluon strength parameter, although the individual contributions to the scattering length have 
important uncertainties so this parameter is not very well determined. In our subsequent 
analysis of phase shifts we shall find that a somewhat smaller value of / ~ —2.5 appears 
reasonable over a wide range of energies for tctt and Kit elastic scattering. 

4. Comparison to Experimental Phase Shifts 

a) I = 7T7T and I = 1/2 Kn S-waves 

We begin our detailed comparison with experiment with a study of the / = vrvr and 
/ = 1/2 Ktt S-waves. We combine these because we expect them to be closely related by 
SU(6) symmetry, and together they provide data on nearly elastic S-wave Ps-Ps scattering 
amplitudes from invariant masses of ~ 0.3 to 1.6 GeV. 

We impose unitarity on the total scattering amplitude by assuming that the phase shifts 
for the individual processes add, as they would for small amplitudes; 

Stot = ^Born + 5r . (35) 

This is equivalent to the prescription used by the LASS collaboration, which was to assume 
a relative phase between resonant and background amplitudes, chosen so that the complex 
sum remains on the unitarity circle, 

sin(5toi)e*'^- = [sm{6,kg)e'''''^] + e^''^"^ ■ [sin(5^)e*'^«] . (36) 

For these amplitudes the Born phase shift is given by (30) and the /o(1400) and -K'q (1430) 
resonance phase shifts by (15-17). We fix the quark model parameters at a^ = 0.6, rriq = 0.33 
GeV, p = 0.677 and p.„ = (3k (from our previous study of / = 2 tttt and / = 3/2 
K-K [jl^] phase shifts). We generally set the pion and kaon masses equal to the isospin 



averages m,r = 0.138 GeV and niK = 0.495 GeV, but for the special case / = vrvr we use 
rriTr = ^«7r+ = 0.1396 GeV, as appropriate for the near-threshold points from Ke^ decay. 



The / = TTTT channel is the most interesting historically ||31[, due to long-standing 
uncertainties in the properties of the very broad /o resonance seen in this slowly rising 
phase shift. This channel is further complicated by inelasticities associated with the narrow 
/o(975), and a scalar glueball is expected near 1.5 GeV, so in this single channel study we 
consider only invariant masses below 0.9 GeV. For our J = vrvr data set we use the s- 



channel and t-channel extrapolations of Estabrooks and Martin |^ (38 points from 0.51 
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GeV to 0.89 GeV), the "case 1" phase shift of Protopopescu et al. |Q (17 points from 0.55 
GeV to 0.89 GeV), and the low-energy data from K^i decay of Rosselet et al. [^ (5 points 
for ^^=0.^=0 _ ^^"=1/=! from 0.289 GeV to 0.367 GeV). We added an estimated low-energy 
P-wave phase shift proportional to k^ to the Rosselet data, with a coefficient chosen to give 
9.4° at 0.51 GeV, as reported by Estabrooks and Martin [^ 



For our Kn data set we use the 24 data points of Estabrooks et al. |^ (from 0.73 to 
1.30 GeV, the full set of separated-isospin phases quoted) and the 37 data points of Aston et 
al. |2l[| below 1.6 GeV. This cutoff was chosen because the radial i^Q(1950) should become 



important above this mass. Since the LASS collaboration |21| tabulated only the K n 



'-^+ 



phase shift [^ rather than the separated isospin ones, we used our formalism to calculate 
I = 3/2 S-wave shifts as well to fit this data directly. We will show the results of the fit for 
/ = 1/2 phase shifts, and the experimental errors we show in the figure are actually those 
quoted by LASS for the mixed-isospin phase 0. This gives a total data set of 121 points for 
S-wave Ps-Ps scattering, from an invariant mass of 0.289 GeV in tttt to 1.59 GeV in Kn. 

We fitted this full S-wave data set with the single s-channel annihilation strength param- 
eter / and masses and widths for the two broad resonances /o(1400) and i^Q(1430). Since 
the mass of the /o(1400) is poorly determined due to our 0.9 GeV cutoff in the vrvr data, we 
constrained the masses by Mj^ = Mx* —0.12 GeV, as suggested by the masses of other mem- 
bers of the ^Pj SU(3) flavor multiplet. There is some indication of an /o(975) contribution 
to the phase shift near 0.9 GeV, so we also added a conventional elastic Breit-Wigner phase 
shift with M = 0.974 MeV and 1=0.047 Gev (PDG values, both flxed) to the / = vrvr phase 
shift. The optimum values of the parameters were found to be / = —2.573, Mx* = 1.477 
GeV, Tk^ = 0.261 GeV, M/, = 1.357 GeV (constrained by Mk^), and F/, = 0.405 GeV. The 
experimental and fltted phase shifts are shown in Fig. 2 {Kn) and Fig.3 (vrvr), and except 
for some inaccuracies in describing the shape of the i^o(1430) resonance region these four 
parameters evidently give a reasonable description of both amplitudes. 

The resonance widths found in this flt are especially attractive; Tf^ = 0.405 GeV is 
reasonably consistent with the recent experimental values of 0.255(40) |^, 0.268(70) [p3| , 
0.375(61) [^, and 0.398(26) GeV ||25[. The discrepancy between these moderate widths 



(from production of the /q in 77 and PP) and the very large widths seen in scattering can 
thus be understood as due to the additional nonresonant contributions to scattering, which 
make the /o and Kq appear broader than they actually are. It is also reassuring that the 
relative fltted /o and Kq widths are not far from the naive SU(6) expectation of F/g/Fx* = 2 
(assuming only tttt and Kn modes and neglecting phase space differences). 

The fltted mass for the broad /o, Mf^ = 1.357 GeV, compares well with experimental 
values from recent PP annihilation experiments, 1.335(30) p3, 1.345(12) p5[, and 1.374(38) 



GeV [Q. Our fltted Kq mass of 1.477 GeV is signiflcantly higher than the LASS value of 
1.412 GeV, although an independent analysis of the LASS data by Weinstein and Isgur 
01 using a two-channel Schrodinger formalism found 1.47 GeV, essentially identical to our 
result. The discrepancy appears to be due mainly to the use by LASS of the phase shift 
formula (15) with -y/i replaced by Mr; on fltting our Kn S-wave data set with their form we 
flnd Mk* = 1.432 GeV. Although the most reasonable generalization of the Breit-Wigner 
form to a broad resonance is not well established theoretically, we found that (15) with -y/i 
replaced by Mji gives an unsatisfactory flt to the low energy I = nn data (a shoulder is 
predicted near threshold), and for this reason we have used the PDG form (15). 
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The scattering lengths in this fit and the quark Born diagram (nonresonant) and reso- 
nance contributions are as follows: 

a^=0(7r7r) = 0.41 fni = 0.24 fm(nonres.) + 0.17 fm(res.) , (37) 

alr^^^iKn) = 0.33 fm = 0.21 fm(nonres.) + 0.12 fm(res.) . (38) 

Although Pennington and especially Burkhardt and Lowe [^ prefer a value closer to 0.3 
fm for the 1 = vrvr scattering length [jl6[, it is clear from Fig. 3 that our fit (with 0.41 fm) 



gives a satisfactory description of the existing low-energy J = tttt data, and considerably 
improved low energy measurements will be required to distinguish these values. The larger 
value of about 0.472(8) fm cited by Estabrooks et al. ||30| for the I = 1/2 Kir scattering 



length may give a somewhat improved fit to the low-energy Ktt data (see Fig. 2), but a better 
measurement of S-wave phase shifts near and below 0.9 GeV would be required to confirm 
this larger value. If improvements in the low energy data do confirm a value of 0.3 fm for 
/ = TTTT, this could be accommodated in our formalism by reducing the broad /o width to 
about 300 MeV and changing / to —2. It is difficult however to see how a Ktt scattering 
length of 0.47 fm could be fitted simultaneously without unrealistic parameter changes. For 
this reason we suggest that the error in (34) was underestimated, and that a value near 0.3 
fm will be found in a higher-statistics measurement of low energy / = 1/2 Kit scattering. 

To test the stability of the fitted parameter values we performed fits with / fixed and 
the resonance masses and widths free (but with Mk* — Mf^ = 0.12 GeV constrained). The 
variation of the fit residual F = fi;I]i(^i(expt.) — 5j(thy.))^/ej(expt.)^ and the resonance 
parameters with / is shown in Fig.4. (We normalize F to unity for / = 0.) Evidently 
nonzero / improves the fit considerably, by better than a factor of two relative to / = 0. 
There is a strong correlation between / (which provides a smoothly rising phase shift) and 
the resonance widths; with / = the best fit requires implausibly large widths of Tj^ = 1.47 
GeV and F^* = 0.57 GeV. As the nonresonant scattering is increased with — /, the widths of 
the resonances required to describe the balance of the low-energy scattering fall. By / = —2 
the widths have fallen to the more reasonable values F/^ = 0.61 GeV and Tk* = 0.34 GeV. 
There is a region of comparable quality of fit, — 2 > / > —3, although by / = — 3 the fitted 
/o width is 0.27 GeV, near the lowest experimental estimates. 

Finally, we tested some of the uncertainties in modelling broad resonances discussed in 
Sec. 2. 3 by carrying out fits to the Kn data set alone with different forms for the resonance. 
The PDG form (15) gives {MK*{GeV),TK*{GeV)J) = (1.477,0.257,-2.602); the LASS 
form (as discussed above) gives (1.432,0.308,-1.403); the LASS form with a single time 
ordering (which is a nonrelativistic Breit-Wigner with the energy-dependent width (16)) gives 
(1.457, 0.284, -2.011); the LASS form with no meson form factor gives (1.477, 0.268, -2.584); 
and finally, a simple nonrelativistic Breit-Wigner phase shift with no meson form factor and 
a constant width gives (1.471,0.290, —2.393). Evidently the various assumptions about the 
resonance phase shift lead to variations of about ±20 MeV in the mass and width and about 
±0.5 in the parameter /. 

h) I = 1 TTTT and 1 = 1/2 Kti P-waves 

The higher partial waves may be treated similarly, and there are fewer complications 
because the resonances are narrower and the Born contributions are smaller. Indeed, a 
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casual inspection of the J = 1 vrvr and / = 1/2 Kit P-wave phase shift data to 1.5 GeV shows 
clear evidence for the p(770) and i^*(892) and little else. 

For our 7 = 1 tttt P-wave data set we again use the s-channel and t-channel extrapolations 
of Estabrooks and Martin |^ (42 points from 0.51 GeV to 0.97 GeV) and the "case 1" phase 
shift of Protopopescu et al. [Q (26 points from 0.55 GeV to 1.15 GeV). For Kt^ we choose 
to fit only the Estabrooks et al. |30| P-wave data (24 points, at energies described above); 
since the / = 3/2 P-wave is not well established we cannot use the mixed-isospin P-wave 



phase shift tabulated in the LASS data ^ 



To fit these P-waves we again use the resonant phase shift (15-17) and set £ = 1 in 
(30). The quark model parameters a^, mg, m^ and (3t, = (3k and the meson masses m.„ 
and rriK were assigned the same values as in the S-wave fit. Since the quark Born diagram 
contributions to the P-wave phase shifts are rather small (typically about 5° at 1 GeV 
invariant mass), we do not fit them independently but instead assume the value of / found 
in the S-wave scattering amplitudes, / = —2.573. This leaves the p and K* masses and 
widths as free parameters. Fitting the 92 P-wave phase shift data points gives Mp = 0.7703 
GeV, Tp = 0.1563 GeV, Mk* = 0.8950 GeV and Tk- = 0.0544 GeV; these are in quite good 
agreement with PDG values, the largest discrepancies being about 5 MeV in both widths. 
Note that the nonresonant amplitudes play an important part in bringing about this close 
agreement; imposing / = makes the fit less than half as accurate (it more than doubles 
F) and leads to p parameters of Mp = 0.7633 GeV and Tp = 0.1434 GeV, which do not 
reproduce the PDG values of Mp = 0.7681(5) GeV and Tp = 0.1515(12) GeV so well. 

The general features of the vrvr and Kn fits are quite similar, so we display results from 
K-K only, which has smaller errors at the higher energies. The fit to the / = 1/2 Kn P-wave 
is shown in Fig. 5, which also shows the individual contributions from the i^*(892) and the s- 
and t-channel Born diagrams. Evidently the quark Born diagram contribution is dominated 
by s-channel gluon exchange (hybrid diagrams). The t-channel diagrams give zero for 7 = 1 
TTTT scattering since they are even under 6 — >■ —6. The Kit P-wave is related to the tttt 
P-wave by SU(6) symmetry, and the Kn t-channel contribution is nonzero only because of 
SU(6) violation through m^ 7^ rnu,d- It appears that the resonance shape alone does not 
give a good description of the phase shift somewhat above the K* (note especially the 0.95 
to 1.25 GeV region, and the Born diagrams (with no free parameters) give just about the 
contribution required for good agreement with experiment. This may be correct, but we 
must be cautious in this interpretation; the smaller resonance contribution is due to the 
energy-dependent width T{E) in (15), which is increasing quite rapidly for the K* in this 
mass region. A similar effect is seen in the / = 1 tttt fit. Since this T{E) is rather uncertain, 
we can only conclude that the quark Born diagrams give a contribution to P-wave phase 
shifts which is consistent with data, and may be observable somewhat above the p and K* 
masses. The systematic uncertainties in describing the resonance phase shifts however are 
comparable to the Born diagram contributions, and until the resonance contributions can 
be established with better accuracy we cannot claim to have confirmed the predicted quark 
Born amplitudes in this P-wave data. 
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5. Conclusion and Suggestions for Future Studies 

We have extended the constituent interchange model of Barnes and Swanson to processes 
with valence qq annihilation by incorporating s-channel gluon exchange and s-channel rela- 
tivized Breit-Wigner resonances. We applied these techniques to J = and 1 = 1 irir and 
/ = 1/2 Kit scattering in S- and P-waves, since these amplitudes are well established experi- 
mentally and can be studied for evidence of nonresonant scattering in addition to the known 
s-channel resonances. In a simultaneous fit to the vrvr and Ktt S-waves we determined the 
strength of the s-channel gluon exchange (hybrid) diagram and fitted the masses and widths 
of the if* (1430) and the problematical /o(1400). In our best fit we find an /o mass and 
width of 1357 MeV and 405 MeV, comparable to recent results from 77 and PP production 
experiments. We conclude that most of the low energy S-wave scattering in these channels 
is due to nonresonant s-channel gluon annihilation, which makes the resonances appear very 
broad in phase shift data. We also applied these techniques to the / = 1 tttt and 7 = 1/2 Ktt 
P-wave phase shifts, and concluded that the s-channel resonances p(770) and if* (892) alone 
suffice to describe most of the P-wave phase shifts from threshold to over 1 GeV. There may 
be evidence for the rather small contribution expected from the quark Born diagrams in the 
P-waves near and above 1 GeV. 

We repeatedly found that uncertainties in the generalization of the nonrelativistic Breit- 
Wigner phase shift to broad resonances limited our ability to separate the resonant and 
nonresonant contributions to scattering in these channels. For this reason it would be very 
useful to establish the limitations and range of validity of the PDG form (15) we have 
assumed in most of this work. 

These techniques could be applied most usefully to reactions in which nonresonant s- 
channel gluon exchange dominates s-channel qq resonances. One such possibility is PP -^ 
AA R^, which has been the subject of detailed experimental investigation at LEAR. It would 



be interesting to determine whether the s-channel gluon Hamiltonian (14), with strength / 
fitted to Ps-Ps elastic scattering, also gives an accurate description of this qq -^ ss annihi- 
lation process. This study would presumably require incorporation of initial- and final-state 
interactions in a coupled channel formalism. This rection has already attracted considerable 
theoretical interest, and analyses using similar techniques have been reported |]3j 
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